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Abstract 

Motivated by potential applications in architecture, we study Darboux cyclides. 
These algebraic surfaces of order < 4 are a superset of Dupin cyclides and quadrics, 
and they carry up to six real families of circles. Revisiting the classical approach 
to these surfaces based on the spherical model of 3D Mobius geometry, we provide 
computational tools for the identificiation of circle families on a given cyclide and 
for the direct design of those. In particular, we show that certain triples of circle 
families may be arranged as so-called hexagonal webs, and we provide a complete 
classification of all possible hexagonal webs of circles on Darboux cyclides. 

Key words: architectural geometry, Mobius geometry, geometry of webs, web 
from circles, Darboux cyclide 



1 Introduction 



The development of Computer- Aided Geometric Design has largely been driven 
by applications. Its origins trace back to the need for computationally efficient 
design and processing of freeform surfaces in the automotive, airplane and ship 
industry. Many other applications (industrial, medical, scientific data analysis 
and visualization) shaped the development of the field. However, architecture 
- which has been using CAD since decades — only very recently attracted 



the attention of the Geometric Modeling community (see Pottmann et al. 



(2007a)). The main reason for the growing interest of Geometric Modeling in 
problems from architecture is the trend towards freeform structures which is 
very clearly seen in contemporary architecture and led by star architects such 
as Frank Gehry or Zaha Hadid. 

While digital models of architectural freeform surfaces are easily created using 
standard modeling tools, the actual fabrication and construction of architec- 
tural freeform structures remains a challenge. In order to make a freeform 
design realizable, an optimization process known as rationalization has to be 



Preprint submitted to Elsevier 



27 January 2013 



applied. This means its decomposition into smaller parts, thereby meeting 
two competing objectives: feasibility, and consistency with the designer's in- 
tentions. Most rationalization techniques replace smooth surfaces, possibly 
with an additional curve network (of panel boundaries) on them, by other 
structures like meshes with special properties. The guiding thought in all con- 
siderations is the efficient manufacturing of the surface parts (panels) and their 
respective necessary supporting/connecting elements. Both simple geometry 
and repetition of elements contribute to this goal of efficiency (see e.g. Liu 



et al. 



et al. 



(2006); Pottmannet al. (2007b, 2008); Schiftner et al. (2009); Eigensatz 



(2010); Fu et al. (2010); 



Singh and Schaefer (2010)). 



As a contribution towards rationalization of architectural freeform structures, 



Bo et al. (2011) have recently studied so-called circular arc structures. A cir- 



cular arc structure (CAS) is a mesh whose edges are realized as circular arcs 
instead of straight line segments and which possesses congruent nodes with 
well-defined tangent planes. Figure [T] shows two examples of circular arc struc- 
tures with quad mesh and triangle mesh combinatorics, respectively. The nodes 
are congruent and possess edge angles of 90 and 60 degrees, respectively. 




Fig. 1. Circular arc structures according to Bo et al. (2011 ). Left: Quadrilateral CAS 



with an edge angle of 90 degrees covering an architectural design. Right: Triangular 
CAS (edge angle 60 degrees) on the Eindhoven Blob by M. Fuksas. 

A circular arc structure (CAS) of regular quad mesh combinatorics consists of 
two discrete sets of arc splines (tangent continuous curves formed by circular 
arcs) which intersect each other under constant, in particular right angle. For a 
triangular CAS, we have three sets of arc splines meeting in regular 60 degree 
nodes. Three families of curves which are the image of a regular triangular grid 



under a diffeomorphism form a so-called hexagonal web or 3-web (see Blaschke 



and Bol (1938)). Thus, a triangular CAS is a special discrete 3-web. 



Now we leave aside the congruence of nodes, but ask for CAS whose generating 
curves are single arcs and not just arc splines. This amounts to the following 
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basic questions: 



(1) Which structures can be built from two families of circles (but not nec- 
essarily meeting under a constant angle) and how can we design them? 

(2) How can we characterize and design 3-webs formed by circles (C-webs)? 

It is conjectured that a surface which carries three families of circles is a 
Darboux cy elide (defined in Section [2j; a proof is missing so far. However, it 
seems to be reasonable to assume that the conjecture is correct; there is no 
known surface different from a Darboux cyclide which carries three families 
of circles. There are, however, surfaces with two families of circles which are 
not Darboux cyclides: an example is furnished by a translational surface with 
circles as profile curves. To keep the present paper more focused and since 
we are mostly interested in webs from circles, we will only study Darboux 
cyclides. Knowing that a complete answer to the first question would require 
additional efforts, we may return to this topic in another publication. 

In the present paper, we will lay the geometric fundamentals for an under- 
standing of Darboux cyclides (which hopefully will also serve as an introduc- 
tion to basic concepts of 3D Mobius geometry) and then proceed towards 
design methods for Darboux cyclides, their families of circles and C-webs. 
Our main result is a classification of all possible C-webs on Darboux cyclides 
(Theorem 17). 



1.1 Previous work 



Darboux cyclides are a topic of classical geometry. Main contributions are 



due to Kummer (1865) and Darboux (1880). Probably the most complete 



discussion is found in the monograph by Coolidge (1916). After a long period 



of silence on these surfaces, geometers again got fascinated by them, especially 



since they can carry up to six families of real circles (see Blum (1980); Takeuchi 



(2000)). Related studies are those which show under certain assumptions that 



a surface with a higher number of circle families must be a sphere (see Takeuchi 



(1985, 1987)). 



In geometric modeling, Dupin cyclides and various extensions attracted a lot 
of attention (see the survey by Degen (2002)), but to our best knowledge, 
Darboux cyclides have not yet been studied from the geometric design per- 
spective, with the noteable exception of a very recent paper by Krasauskas 



and Zube (2011), whose relation to cyclides has been realized in a discussion 



of our research with Rimas Krasauskas. 



The problem of determining all hexagonal webs from circles in the plane (or 
equivalently, on the sphere) turned out to be very difficult. A special class 
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of webs from circles, initially found by Volk (1929) by analytical methods 



has soon been realized as a simple projection of webs from straight lines in 



the plane onto the sphere (Strubecker (1932)). An elegant construction of a 



nontrivial class of webs from circles by Wunderlich ( 1938 ) will turn out to be 



an important ingredient in our study of C-webs on Darboux cyclides. Recently, 



Shelekhov (2007) could classify all webs formed by pencils of circles. It is still 



unknown whether the so far described instances of webs from circles are the 
only possible ones. 

Overview. Our paper is organized as follows. Section [2] studies Darboux cy- 
clides by employing the spherical model of Mobius geometry and mostly meth- 
ods of projective geometry (pencils of quadrics) to get insight into the circle 
families they contain. Based on these insights, we are able to devise design 
methods for Darboux cyclides (Section [3]) and to characterize all C-webs on 
them in Section |4} Finally, we address future research directions and open 
problems in Section [5] 



2 Darboux cyclides and their families of circles 



We start this section with the definition of a Darboux cyclide. In Section 2.1 we 



introduce the spherical model of Mobius geometry and in Section 2.2 we show 
that a Darboux cyclide is represented by the carrier of a pencil of quadrics 
in this model. Using this representation, we describe families of circles on a 



Darboux cyclide in Section 2.3 and determine the possible number of such 



families in Section 2.4 In Section 2.5 we apply polarity in the spherical model 



to set up some background for design of Darboux cyclides in the next section. 

A Darboux cyclide is a surface whose equation in a Cartesian coordinate system 
has the form 

D : \(x 2 + y 2 + z 2 ) 2 + (x 2 + y 2 + z 2 )L(x, y, z) + Q(x, y, z) = 0, (1) 

with a constant A, a polynomial L = fix + vy + kz, and a polynomial Q of 
degree at most 2, not vanishing simultaneously. Hereafter the coefficients and 
the coordinates are real numbers unless otherwise explicitly indicated. 

If the left-hand side of equation ([I]) factors into non-constant polynomials in 
x, y, z with (possibly complex) coefficients, then the Darboux cyclide is called 
reducible. A reducible cyclide either splits into a union of spheres/planes or 
degenerates to a curve in IR 3 . A cyclide of the form a(x 2 + y 2 + z 2 ) + bx + cy + 
dz + e = is itself a sphere/plane and is referred as trivial hereafter. Thus, in 
the following we focus only on irreducible and nontrivial cyclides. 

The intersection curve of a cyclide D with a sphere S : x 2 = Li(x), where 
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x := (x, y, z) and L\ is a polynomial of degree at most 1, lies on the quadric 

Qf. AL 1 (x) 2 + L 1 (x)L(x) + Q(x) = 0, 

and thus it is in general an algebraic curve c of order at most 4, which is 
common to all quadrics of the pencil spanned by S and Q\. Hence, if the 
sphere S touches D at two points, the curve c has two double points and 
thus it is formed by two circles (or a doubly counted circle c, in which case 
the sphere touches D along c). This already indicates the presence of circles 
on cyclides, which is better studied by a slightly more advanced method as 
follows. 

2. 1 Some basic facts on the spherical model of Mobius geometry 

We map points x = (x,y,z) G 1R 3 to points x = (x\, . . . ,X4) G M 4 on the 
sphere £ : x\ + . . . + x\ = 1 via the inverse stereographic projection. The 
projection center shall be Z = (0, 0, 0, 1) and we embed 3-space into 4-space 
via x i — y (x,y,z,0) G M 4 . Corresponding points x and x lie on straight lines 
through Z. The mappings between x and x in both directions are described 
by the formulae 

x = (xi, . . . , x A ) = 1 (2s, 2y, 2z, x 2 - 1), (2) 
x^ + 1 

x = (x, y, z) = — - — (x 1 ,x 2 ,x 3 ). (3) 

1 — X4 

By definition, the set of points in Mobius geometry is Euclidean 3-space with 
one ideal point oo added; the ideal point oo corresponds to the projection 
center Z. 

Unfortunately, it is hard to illustrate geometric facts involving a 4D space. 
Thus, to support the reader's imagination, Fig. [2] illustrates the lower dimen- 
sional counterpart of the stereographic projection of a sphere in 3-space onto 
a plane. What we are going to say below about spheres, is true for circles in 
the lower dimensional counterpart. 

A Mobius sphere (or M-sphere) S is the set (strictly speaking, in a complex 
3-space) given by the equation Ax 2 + ax + by + cz + d = 0, where A, a, b, c, d do 
not vanish simultaneously. By definition, oo G S if and only if A = 0. A Mobius 
sphere is either a plane (A = 0), or a Euclidean sphere (real sphere henceforth), 
or a point (null sphere), or does not have real points at all (imaginary sphere). 

One can work with imaginary spheres without appealing to complex numbers 
as follows. The equation of a nonplanar M-sphere can be rewritten in the form 
(x — m) 2 = r 2 , where m G R 3 is the midpoint and r 2 G IR is the square of the 
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Fig. 2. The spherical model of planar Mobius geometry and stereographic projection: 
any planar section S of the sphere £ corresponds to an M-circle S and the extended 
stereographic projection maps the pole S* of the plane to the center m of the circle. 
The projection center Z £ £ corresponds to the ideal point oo and circles through 
Z correspond to straight lines in the plane. 

radius. In the case of imaginary sphere S the number r 2 is negative; we can 
write r 2 = — p 2 for some p > 0. Define the real representer S r of the sphere 
S as the sphere with center m and radius p. In the following, all statements 
involving imaginary spheres can be easily converted to statements on their 
real representers using Remark [T] below. 

M-spheres correspond to intersections S of E with hyperplanes 

H s : ax\ + bx 2 + cx 3 + (A — d)x± + A + d = 0; 

of course, the 2-surfaces S are M-spheres, possibly imaginary or null. If the M- 
sphere S is a plane, it contains the ideal point, i.e., H s and S pass through Z. 
Often one uses the polarity with respect to E. Using homogeneous coordinates 
X := (Xl, . . . ,X 5 ) in the projective extension P 4 of the 4-space, with Xi = 
Xi/X 5 , i = 1, . . . , 4, the pole S* of Hs is given by 

S* = (a, b, c, A — d, —A — d). 

The coordinates of S* are called pentaspherical coordinates of the M-sphere S. 
For an M-sphere S with center m = (m x ,m y ,m z ) and squared radius r 2 , we 
have 

A = 1, a = —2m x , b = —2m y , c = —2m z , d = m 2 — r 2 . (4) 

The stereographic projection is strictly speaking a mapping from points of the 
sphere E C M 4 to 1R 3 U {oo}. However, it is useful to also use the underlying 
central projection with center Z and image space X4 = for our considerations, 
even if we do not apply it only to points of E; we will call this projection the 
extended stereographic projection henceforth. From Q we can easily verify the 
well-known fact that the extended stereographic projection maps the pole S* 
onto the midpoint of S (see also Fig. [2]). 
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An M-circle k of 3D Mobius geometry is the intersection of two M-spheres 
(strictly speaking, in the complex 3-space) and thus it is either a circle, a 
straight line, a single point, or the empty set from a Euclidean perspective. 
Formal definitions of real, null, and imaginary circles are left to the reader. 
An M-circle is mapped to the intersection k of E with a 2D plane, i.e., to an 
M-circle in E. The polarity with respect to E maps the plane of k to a straight 
line k*. The points S* of k* determine M-spheres S through k. These spheres 
are said to form a pencil of spheres. 

A pencil is more generally defined as the set of all M-spheres S whose cor- 
responding points S* form a straight line k* in P 4 . If k* does not intersect 
E, we obtain an elliptic pencil, formed by all M-spheres through a real cir- 
cle. If k* is tangent to E, we get a parabolic pencil, and if k* intersects E in 
two points, we have a hyperbolic pencil. A parabolic pencil is formed by all 
M-spheres touching a plane P at a fixed point p G P. The common set of a 
hyperbolic pencil is an imaginary circle; the pencil contains two null spheres, 
corresponding to the points E fl k*. Clearly, in all cases the centers of the M- 
spheres lie on a straight line A, which is also the image of the line k* under the 
extended stereographic projection. Due to the rotational symmetry of pencils, 
it is sufficient to understand the planar sections through the rotational axis 
A, which are pencils of circles (see Fig. [3]). 



Fig. 3. Pencils of circles in 2D Mobius geometry are stereographic projections of 
pencils of circles 5 on E; the latter are easy to understand since they are cut out 
from E by a pencil of planes. The poles S* of these planes lie on a straight line k*. 
The figure illustrates the 3 types of pencils, which by rotation about the central line 
generate the 3 types of pencils of spheres. 

Linear families of spheres S are those where the set of poles S* is a d- 
dimensional subspace of P 4 . Apart from the already discussed pencils (d = 1), 
we still have to consider the cases d = 2 (bundles of spheres) and d = 3 (sphere 
complexes). 
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Before doing so, we should look at orthogonal M-spheres, i.e., those which 
intersect under a right angle. Two real spheres S\, S2 with centers mi, m 2 and 
radii r 1; r 2 , are orthogonal if (mi — m 2 ) 2 = t\ + r|. Their corresponding points 
S* (in pentaspherical coordinates) are 

S* = (2m ijX ., 2m i)2/ , 2m itZ , m 2 - r 2 - 1, m 2 - rf + 1). 

Two points A = (ai, . . . , a 5 ) and 5 = . . . , 65) in P 4 are conjugate with 
respect to S (one lies in the polar hyperplane of the other one) if (A,B) := 
di&i + . . . + 0,464 — a 5 6 5 = 0. We see immediately that (S^, S 2 ) = if and only 
if (mi — m 2 ) 2 = r\ + r\. Hence, orthogonal real spheres Si, S* 2 are character- 
ized by conjugate points S 1 *,^. The latter characterization allows to define 
orthogonality also for null and imaginary spheres. 

Remark 1 In the above discussion, it can happen that one of the spheres, say 
S 2 , is imaginary, and thus corresponds to a point S 2 G M 4 in the interior of the 
3-sphere S. Hence, the orthogonality condition now reads (m 1 -m 2 ) 2 = r\ — p\, 
where p\ := —t\. Then, the two spheres S\ and S2 are orthogonal if and only 
if the intersection circle of Si and S 2 is a great circle of S 2 (see Fig. 1^1). 




Fig. 4. Orthogonality between M-spheres S±,S2 is not only defined for real spheres 
(left), but also for imaginary ones. A real sphere Si is orthogonal to an imaginary 
sphere S* 2 (with center m 2 , radius ip-z), if it intersects the real representer S 2 (with 
center m 2 , radius p 2 ) of S2 in a great circle (right). 

It is now easy to understand a linear complex of spheres. The points S* cor- 
responding to the M-spheres of the complex lie in a hyperplane H, whose 
pole shall be denoted by A*. The point A* is conjugate to all points S* C H 
and it defines a (not necessarily real) M-sphere A, which is orthogonal to all 
M-spheres S of the complex (and which is seen in the spherical model as inter- 
section A = H n S). Hence, a complex is formed by all M-spheres orthogonal 
to a fixed M-sphere A. We have again 3 types, depending on whether A is real, 
imaginary or null, i.e., degenerates to a point. In the latter case, all M-spheres 
of the complex pass through this point. 

A bundle of spheres is defined by a 2D plane P C P 4 , which may be viewed as 
the intersection of two hyperplanes Hi, if 2 . Hence, the M-spheres of a bundle 
belong to two complexes and thus they are orthogonal to two (not necessar- 
ily real) M-spheres Ai,A 2 . This requires orthogonality to the (not necessarily 
real) intersection M-circle A\ D A 2 , and thus the M-spheres of a bundle have 
their centers in a plane and they are all orthogonal to an M-circle k in that 
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plane (note that the M-circle k is seen in the spherical model as the inter- 
section M-circle k = PflS). Once again there are three types, since k may 
be real, imaginary or null (see Fig. [5J. The case of a imaginary circle is easy 
to understand, since now all M-spheres of the bundle pass through two real 
points (on the real axis of the imaginary circle). This follows for example from 
the fact that in this case two hyperplanes through P are tangent to E (or by 
observing that the polar dual to the plane P is a straight line, determining 
a hyperbolic pencil of spheres which are all orthogonal to the spheres of the 
bundle). 



Fig. 5. The three types of sphere bundles are defined as sets of M-spheres which 
intersect an M-circle k under right angle. Top left: for an elliptic bundle, k is real. 
Bottom: in a hyperbolic bundle, k is imaginary, implying that the spheres pass 
through 2 fixed points. Top right: a parabolic bundle is a limit case and formed by 
all M-spheres centered in a plane and passing through a fixed point in that plane. 

Mobius transformations (or M-transformations) act bijectively on the points 
and spheres of Mobius geometry. M-transformations are the only conformal 
(angle preserving) maps in Euclidean 3-space; they include Euclidean similar- 
ities as special cases. In the spherical model they are seen as those projective 
transformations k: P 4 — > P 4 which map £ onto itself. Using homogeneous 
coordinates, the equation of £ is X T • J ■ X = 0, with the diagonal matrix 
J = diag(l, 1, 1, 1, —1). A projective map k : X i— y A ■ X in 4-space maps S 
onto itself and thus induces an M-transformation if and only if its matrix A 
satisfies 



for some A ^ 0. 

The basic example of an M-transformation is the inversion with respect to an 
M-sphere S. In the spherical model it appears as a projective symmetry of 
E, more precisely, a harmonic perspective collineation k which maps S onto 
itself. The axis of this collineation is the hyperplane H s D S and the center is 
its pole S*. Thus an inversion with respect to an imaginary sphere 5 is also 




A T ■ J ■ A = A J, 



(5) 
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defined. In this case, S* lies in the interior of E. In 3-space, such an inversion 
can be realized by reflection at the (real) center of S and subsequent inversion 
with respect to the real representer S r of S. Notice that the real representer 
is not M-invariant, i.e., K(S r ) ^ (^S) r in general. 



2.2 Darboux cyclides in the spherical model of Mobius geometry. 



In this subsection we will elaborate on the classical result that a Darboux 
cyclide D is the stereographic projection of a surface D in S which is obtained 



as the intersection o/E with another quadric T (see Coolidge (1916); here by 
definition the ideal point oo belongs to a nontrivial cyclide D if and only if 
A = in equation ([!])). In other words, D is the common set of points, the 
so-called carrier, of a pencil of quadrics which contains the sphere E. Deriving 
this result is not hard, and it will give us important insight needed later in 
our studies. 

Let us first show that inverse stereographic projection maps a Darboux cyclide 
D to the carrier of a pencil of quadrics. For that, we substitute from ^ into 
(|T|. Splitting Q into homogeneous terms of degree 2, 1, and as Q(x) = 
g(x) + q ■ x + g and writing L — 1 • x yields with x := (xi, x 2 , £3) 

/ x 2 \ 2 x 2 l-x q(x) qx 

+ T1 ^ + n \? + ^ — + qo = o. 



\{1 — X4) 2 J (1 — X4) 2 1 — £4 (1 — X4) 2 1 — £4 

Inserting the equation of E, x 2 = 1 — x 2 , we finally obtain the equation of a 
quadric or a plane in 4-space, 

T : A(l + x 4 ) 2 + (1 + £4)1 ■ x + q(x) + q • x(l - x 4 ) + q (l - x 4 ) 2 = 0. 

We may assume without loss of generality that T is a quadric by adding a 
multiple of x 2 + x\ — 1 to the left-hand side of the equation. The quadric T 
contains the projection center Z if and only if A = 0, i.e., 00 G D. 

Conversely, we have to show that intersection of E with any quadric T C M 4 
yields a surface D whose stereographic image in 3-space is of the form (jl|, i.e., 
a Darboux cyclide. Again separating terms of different degrees and splitting 
(xi, . . . , x 4 ) as (x, £ 4 ), we write the equation of the quadric as 

r : q 2 (x) + q\{x)x A + q$x\ + h(x) + / ^4 + c = 0, 

where q 2 , {qi, h} and {go, ^0, Q)} are of degrees 2, 1, and 0, respectively. Inser- 
tion from ([2]) yields 

(go + ^0 + c )(x 2 ) 2 + 2x 2 [gx(x) + Z x (x)] + 4g 2 (x) + Q(x) = 0, 
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with 

Q(x) = 2(c - g )x 2 + 2(/i(x) - gi(x)) + q - k + c . 
This equation is of the form and thus represents a Darboux cyclide. 

Transferring the study of Darboux cyclides to the study of the carriers of 
pencils of quadrics which contain the sphere £ has a lot of advantages since 
we are now able to apply knowledge from projective geometry on pencils of 
quadrics. The study of all possible types of cyclides, i.e., pencils of quadrics 
in 4-space, can be elegantly based on the Weierstrass theory of elementary 
divisors (a workaround being found 
p. 301) andp5ij|([l976|). 



1916 



Takeuchi 


( 


2000 


); see also ( 


Coolidge 



2.3 Families of circles on Darboux cyclides. 



Working in the spherical model, we represent the Darboux cyclide as the 
carrier D of the pencil fP of quadrics 

X T ■ (A — tJ) • X = 0, (6) 

spanned by the sphere S : X T - J-X = 0, and another quadric V : X T -A-X = 0. 
We can work with the inhomogeneous pencil parameter t, since we do not have 
to represent E. Obviously, given a point Y ^ £ (with coordinates Y), there is 
a unique quadric Ty in <P which contains Y, namely the one with the pencil 
parameter t = (Y T ■ A ■ Y)/(Y T • J ■ Y). 

If the carrier D of !P contains a circle k in a plane Pk, we can choose any 
point Y ^ k in Pk and investigate the quadric Ty £ & which contains Y. The 
intersection of the plane Pk with Ty must contain the circle k and the point Y 
and thus the entire plane Pk must be contained in the quadric Ty. A regular 
quadric in 4-space can contain at most straight lines, but no planes. Hence, 
Ty is singular, i.e., a quadratic cone. 

Hence, looking for circles of D, we have to search for their planes on the 
quadratic cones in fP. This is a well-known subject, and can be analytically 
accessed as follows. A quadric in the pencil for some parameter t is a cone if 

det{A - tJ) = 0. 

Hence, the possible parameter values ti are the real eigenvalues of the matrix 
A ■ J~ x = A ■ J. Let ti be such an eigenvalue. Then the linear system (A — 
ti J) ■ X = has a nontrivial solution, say Vj, which describes a point Vi, the 
vertex of the cone. This implies 

A-Yi = tiJ ■ Vj, 



11 



which says that the polar hyperplanes of the point Vi with respect to the 
quadric rGf:X'-i-X = and with respect to the sphere £ agree. In 
fact, the polar hyperplanes of V, with respect to all quadrics in 2> are identical, 
because of 

(A-tJ)-Vi = (ti-t)J-Vi. 

We see that Vi is a vertex of the common polar simplex of the pencil fP (which 
needs not be a regular simplex, needs not have only real vertices and may 
degenerate in various ways). In the generic case, we get 5 linearly independent 
eigenvectors of the matrix A ■ J, determining the 5 vertices of the common 
polar simplex. Only real vertices Vi are of interest to us. Connecting any of 
these vertices V, with the carrier D, we obtain a quadratic cone I\ which 
contains D. 




Fig. 6. A cyclide D determines a pencil of quadrics in 4-space. Each vertex Vi of 
the common polar simplex of the pencil corresponds to a sphere Bi in 3-space with 
respect to which D is symmetric (inversion with respect to Bi maps D onto itself). 
This figure shows a cyclide D which is symmetric with respect to 5 (pairwise orthog- 
onal) spheres B\ , . . . , B5 ; since B$ is imaginary, we illustrate its real representer B§ . 

Remark 2 Any real vertex Vi ^ S of the common polar simplex, together with 
its polar hyperplane, define a joint projective symmetry of all quadrics in the 
pencil, and thus of its carrier D. The vertex Vi determines an M-sphere, say 
Bi, in 3-space (in the notation used above, we have B* = Vi) . The inversion 
with respect to such a sphere Bi maps the cyclide D onto itself. Since the 
vertices Vi are pairwise conjugate with respect to E ; the spheres Bi are pairwise 
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orthogonal. In the case with 5 real vertices Vi, we get cyclides D which are M- 
symmetric with respect to 5 pairwise orthogonal M-spheres (see Fig. [#| and 
Coolidge (1916)); note that one of the 5 points Vi must lie in the interior of 
E and thus the corresponding sphere is imaginary (B 5 in Fig. 




Fig. 7. A pencil of quadrics in 3-space has an algebraic curve c of order 4 as its 
carrier. There are four cones Tj in this pencil; projecting c from a cone vertex Vi 
onto a plane Hi results in a conic Di. 

Once again we look at a lower-dimensional counterpart to visualize the situ- 
ation. Fig. [7] shows a pencil of quadrics in 3-space, spanned by a sphere and 
another quadric. The pencil contains four cones with the vertices at the ver- 
tices of the common polar simplex. The figure also shows that the projection 
of the intersection curve (algebraic curve c of order four) from a point Vi onto 
a plane is a conic Di and that not all points of Di need to be projections of 
real points of c. 

Let us return to 4-space. We may project D from a cone vertex Vi onto a 
hyperplane Hi (not passing through V) arid obtain a quadric Di C Hi. The 
structure of this projection becomes clear in C 4 . The quadric Di appears as 
"doubly covered", i.e., each of its points is the image of two points of D. 
However, not any real point of the quadric Di is the image of two real points 
on D (see Fig. it can as well be the image of two conjugate complex points 
of D. The double coverage of Di explains why it can be the image of a 2- 
surface D of algebraic order 4 under a central projection. Returning to M 4 
we get that if the cone Tj contains planes, they pass through circles on D 
and of course through Vi and thus these planes (and the circles) are projected 
onto straight lines of the quadric Di. This shows us how to obtain circles on 
Darboux cyclides (see Fig. [8]). We summarize our findings in the following 
theorem which is a classical result, but probably not so straightforward to 



used in later studies on cyclides (Blum (1980); Takeuchi (2000)). 



extract from the classical texts (e.g., Coolidge (1916)) so that it did not get 
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Theorem 3 In the spherical model of Mobius geometry, a Darboux cyclide D 
appears as the carrier D of a pencil fP of quadrics which contains the sphere 
£. Any real vertex of the (possibly degenerate) common polar simplex of the 
quadrics in 'S is the vertex of a quadratic cone through D. If such a cone 
contains planes, these planes intersect £ in circles of D, which determine 
circles or lines on D. Any circle or line on D can be found in this way. 



Fig. 8. Extended stereographic projection of the situation in 4D. (Top left) Con- 
necting the image point V( of the cone vertex Vi with the rulings of the projected 
ruled quadric D[ results in planes T which carry circles of the cyclide D. (Top right) 
Each plane T carries two circles and is tangent to the cyclide D at < 2 points, the 
common points of these two circles. (Bottom left) Each plane T is tangent to a 
quadratic cone Cj with vertex V- which is tangent to D[. (Bottom right) A paired 
family of circles arising from the two families of rulings in D\. 

This result gives rise to several conclusions on circle families on nontrivial 
irreducible Darboux cyclides (containing no spheres or planes as a whole). 
There are two types of quadratic cones in the 4-space which contain planes 
(but no hyperplane): 

(1) The cone I\ is obtained by connecting a nonsingular ruled quadric Di in 
a hyperplane Hi with a point Vi outside Hi. The ruled quadric Di carries 
two families of straight lines {rulings). Each plane through a ruling of 
one family intersects Di in a ruling of the second family. Connecting 
with the vertex Vi and interpreting the scene from the Mobius geometric 
perspective, we find two circle families on D (and D). We call them paired 
families, since one family determines the other one. Any sphere through 
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a circle of one family intersects the cyclide D in a second circle which 
belongs to the other family. This holds in particular for the plane of 
such a circle: it intersects the cyclide in a circle of the other family (see 
Fig. [8]). We shall say that paired families are special, if Vi G £ — this is 



only possible if D is the image of a quadric under an inversion ( Coolidge 
1916, cf. Theorem 35 in p. 301). Paired families were called residual in 



Coolidge (1916). 




Fig. 9. Single families of circles are found only on those Darboux cyclides which are 
envelopes of a one-parameter family of spheres. The spheres touch the cyclide D 
along the circles of the single family, and the planes of the circles pass through a 
fixed line (image of the vertex line Vi of a quadratic cone containing D). 



(2) 



The cone Tj is the connection of an irreducible quadratic cone Di (with 
vertex Wi) in a hyperplane Hi and a point Vi outside Hi. This is the 
same as connecting a conic (on Di) in some plane P C Hi with a straight 
line Vi = ViWi which does not intersect P. The straight line is the vertex 
line of the cone. Such a cone carries a single family of planes (through 
and thus we speak of a single family of circles on D and on the 
underlying Darboux cyclide D. Obviously, a sphere through a circle of 
a single family intersects the cyclide in a circle of the same family, just 
like a plane through a ruling of the quadratic cone Di intersects the cone 
in another ruling. However, there are also tangent planes of the cone Dj. 
They lead to spheres touching the cyclide along the circles of a single 
family and thus we see that a Darboux cyclide which contains a single 
family is a canal surface, i.e., the envelope of a one-parameter family of 
spheres (see Fig. [9|. In fact, almost each such Darboux cyclide is the 



image of a complex quadratic cone under a complex inversion (Coolidge 



1916 Theorem 38 in p. 302). 



Proposition 4 There are two different types of circle families on nontrivial 
irreducible Darboux cyclides D. (i) Paired families are two families such that 
a sphere through a circle of one family intersects the cyclide in another circle, 
which belongs to the second family. The planes of both circle families are ei- 
ther the tangent planes of a quadratic cone or the planes of two pencils with 
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intersecting axes (in the sense of projective geometry), (ii) The circles of a 
single family are the characteristic circles of a generation of the cyclide as a 
canal surface; a sphere through a circle of a single family either touches the 
cyclide along that circle or it contains another circle of D belonging to the 
same family. The planes of the circles lie in a pencil. 

The statements on the planes of circles follow from the extended stereographic 
projection of the plane family in the quadratic cone (see Fig. [8]). Indeed, 
we may assume that Z ^ Vi and Z e Hi. If Z ^ Di in type (i) then the lines 
in the irreducible quadric Di are projected to tangents to a certain conic in 
3-space. Joining the resulting lines with the projection of the point Vi we get 
tangent planes of a quadratic cone. The case of two pencils with intersecting 
axes in type (i) arises if the projection center Z, the vertex Vi, and a point of 
the ruled base quadric D { are collinear. This is only possible if A = in (JIJ); 
an example is the paired family of circles in an ellipsoid (not of revolution). 

Using this approach one can obtain further information about the arrangement 
of circle families on Darboux cyclides. Once the quadric Di contains at least 
one line, it contains a line through each point. Thus any family of circles on an 
irreducible cyclide D covers the entire cyclide. Now assume that there are two 
circle families on D. Then through a generic point p e D one can draw circles 
from both families. There are two possibilities. First, these circles may have 
another common point (or touch each other). In this case they lie in one sphere. 
Hence by Proposition [4] the circles belong to two paired families. Second, the 
point p may be the unique transversal intersection point of the circles. In this 
case any two other circles from the families have also a single intersection 
point because the parity of the number of intersection points (counted with 
multiplicities) does not depend on the choice of two particular circles from the 
families. 

Proposition 5 In nontrivial irreducible Darboux cyclides, two circles from 
distinct families, which are not paired together, intersect each other transver- 
sally at a single point. 

Remark 6 Similarly, there are two types of families of conies in the inter- 
section D of two arbitrary nonsingular 3 -dimensional quadrics in 4-space (D 
containing no 2-dimensional quadrics as a whole). Paired families are such 
that a hyperplane through a conic of one family intersects D at another conic, 
which belongs to second family. A hyperplane through a conic of a single fam- 
ily either touches D or intersects it in another conic of the same family. This 
is proved analogously to Proposition [^} 
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2.4 The number of families of real circles on a Darboux cyclide. 



Let us further discuss the circle families on cyclides and start with some simple 
examples. Darboux cyclides which contain two single families must be the 
envelopes of two families of spheres and thus these are the Dupin cyclides 
which have already found a lot of attention in Geometric Modeling (see, e.g., 
Degen (2002); Krasauskas and Maurer (2000)). As a special case, consider a 
ring torus, which obviously contains two single families (meridian circles and 
the circular paths of the generating rotation); furthermore it carries a paired 
family, namely the two families of Villarceau circles. Likewise any Dupin ring 
cyclide carries four families of circles, since it is the image of a ring torus 
under an inversion. Rotating a circle about an axis (not meeting the axis of 
the circle and not contained in the plane of the circle) we may obtain a ring 
shape, which is a Darboux cyclide carrying five real families of circles (two 
paired ones and the single family of rotational circles; for in illustration, see 



Pottmann et al. (2007a)). 



However, there are also Darboux cyclides which carry six families of real circles. 



This has first been pointed out by Blum (1980), at hand of a special case of 



symmetric cyclides, and in a quite involved elementary proof which gives less 
insight than our approach. Blum mentions that Darboux had pointed to 10 
families of circles, but did not differentiate between real and complex ones. 



Takeuchi (2000) provides a discussion of the real families of circles on cyclides, 



similar to our approach, but he does not address the geometric relation to the 
cones in the pencil of quadrics. In none of these papers, we find a figure of such 
a surface, which of course is easy to generate with our approach and provided, 
e.g., in Fig. [10 



We can easily obtain these results and more as follows. Without caring about 
reality and looking only at the generic case, we get Darboux's 10 families, two 
families from each of the 5 quadratic cones. To determine the number of real 
circle families on a cyclide we have to investigate the quadratic cones in a pencil 
of quadrics through E. Using the notation from equation ([6]), the matrices 
of the quadratic cones Tj are C{ := A — £jj, where £, is a real eigenvalue 
of the matrix A ■ J. A quadratic cone with two families of planes (which 
determines a paired family of circles) is given by a matrix C{ of signature 
(0, +, +, — , — ). For a cone determining a single family, the matrix has signature 
either (0, 0, +, +, -) or (0, 0, +, -, -). 

For brevity, let us now focus on the generic case with 5 distinct vertices of 
the common polar simplex, where we also assume that all five vertices are 
real. In this case, we can apply a projective map in 4-space which maps these 
vertices to the origin and the ideal points of the four coordinate axes. In 
other words, we have achieved that matrix A is a diagonal matrix, say with 
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Fig. 10. A cyclide with 6 families of circles, shown in different colors. We will see 



later that we can form eight different 3- webs from these circles (cf. Figure 18). 

diagonal (a±, . . . , a^). Matrix A ■ J is then also diagonal with eigenvalues U ■ 
{at, . . . , a 4 , — a 5 }, and so the eigenvalues of Cj are (ai — U, . . . , a 4 — t i} a 5 + 



Example 7 To give an example of a cyclide with 6 real circle families (Fig.\lu\), 
we consider the pencil with the matrix A = diag(— 2, — 1, 1, 2, 0), which defines 
a quadratic cone r\ with vertex at the origin; Ti contains two plane families. 
Using Cartesian coordinates, our pencil is spanned by 

xj + ...+x 2 A -l = 0, -2x1 ~ x l + x l + 2x l = °- 

Following the above procedure, or equivalently, eliminating one of the vari- 
ables, we arrive at the following 4 additional cones (cylinders parallel to the 
coordinate axes): 

T 2 (t t = -2) : x 2 2 + 3x 2 3 + Ax l ~ 2 = 0, T 3 (U = -1) : -x\ + 2x\ + 3x\ - 1 = 0, 

r 4 (t< = 1) : -3x1 - 2x\ + x l + 1 = 0. r s(^ = 2 ) : ~^ x \ ~ ^> x l - x l + 2 = °- 
We see that T 3 and T 4 determine paired families, while the largest and smallest 
U determine cylinders with an oval base quadric and thus do not contribute 
further real circles. 

Remark 8 For the smallest real eigenvalue t of the matrix A ■ J the matrix 
A — t J cannot have signature (0, +, +, — , — ). Indeed, the eigenvalues of the 
matrix A — tJ change continuously with decreasing t, but for t < t they 
should have the same signs (+, +, +, +, — ) as those of J. Similarly, the largest 
real eigenvalue does not lead to a paired family. Thus there are at most 3 
paired families . Single families may only appear for eigenvalues of multiplicity 
2, hence we never have more than 6 families of real circles on a nontrivial 
irreducible Darboux cyclide. 
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This can easily be cast into a simple algorithm for computing the circle fam- 
ilies on a given cyclide. Assuming that we have already identified the cor- 
responding pencil of quadrics (JgJ) , we compute the real eigenvalues of the 
matrix A ■ J. For each such eigenvalue U, we compute the eigenvalues of the 
matrix A — tiJ. If the eigenvalues of the matrix A — tiJ have signs either 
(+, +, — , — , 0) or (+, +, — , 0, 0) or (+, — , — , 0, 0) then we bring the matrix to 
the form diag(A^, k 2 , — fcf , — fc|, 0) by a linear transformation k of coordinates, 
where ki,k 3 > and k 2 ,k 4: > 0. Denote by Pi(s) and P 2 (s), where s is a 
parameter, the images of the planes 



under the transformation k -1 . In particular, the families Pi(s) and ^(s) co- 
incide, if either k 2 = or k± = 0. Stereographic projections of Pi(s) fl S and 
P2{s) H S are families of circles on the cyclide. If the eigenvalues of A — t{ J 
have signs neither (+, +, — , — , 0) nor (+, +, — , 0, 0) nor (+, — , — , 0, 0) then t{ 
does not contribute to circle families. 



2.5 The polar dual. 

While not absolutely necessary to get further insight, we think that it is helpful 
to also consider the polarity with respect to X and view the situation from this 
perspective. Strictly speaking, the polarity maps points to hyperplanes, but 
it defines a duality which maps spaces of dimension d to those of dimension 
3 — d for d = 0, . . . , 3; this is why we speak of the "polar dual". An M-sphere 
S gives rise to a point S*, the pole of the hyperplane carrying the sphere 
S C S. An M-circle k in IR 3 determines a circle k C S in some plane, whose 
polar image is a straight line which we denote by k*. Recall that the extended 
stereographic projection maps the point S* which determines a sphere S onto 
the center of S. Likewise, it maps the straight line k* determining a circle k 
to the rotational axis of k. 

A circle family of a Darboux cyclide is completely determined by a quadratic 
cone Ti which carries a family of planes. We now look at the polar dual V* of 
the cone, and as above distinguish the two cases. 

Paired circle families. If Tj contains two families of planes and thus de- 
termines two paired circle families on a cyclide, its polar image V* is a ruled 
quadric in some hyperplane V* (the polar hyperplane of the cone vertex Vj). 
Each ruling of T* is polar to a plane in Tj and thus determines a circle on the 
cyclide. The extended stereographic projection maps the ruled quadric T* to 
a quadric (or a plane) cK 3 . The axes of the circles in the paired families 




and 



k\X\ — sk 2 x 2 — k 3 x 3 — sk^Xi = 0, 
skiXi + k 2 x 2 + sk 3 x 3 — k 4 X4 = 0; 
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are projected to two families of lines in Gi. 



The cone vertex V% represents an M-sphere; let us call it Bi (note: B* = Vi) and 
consider the associated complex of spheres S which are orthogonal to Bi. The 
poles S* representing these spheres S lie in the hyperplane V*. This implies 
(recall Remark [2] on M-symmetries) that the cyclide is symmetric with respect 
to Bi. Note that each of the circles in the two families is also symmetric with 
respect to Bi, i.e., orthogonal to Bi, since it is an intersection of two spheres 
of the complex. 



\ 




Fig. 11. A paired family of circles on a cyclide D can be constructed as follows: 
Prescribe a ruled quadric G and a sphere B, and consider all circles which are 
orthogonal to B and have a rotational axis on G. This may be viewed as a trans- 
formation which maps the rulings of G to circles on D. 

Proposition 9 All circles of two paired families are orthogonal to a fixed Ad- 
sphere Bi and their axes form the two families of rulings on a ruled quadric 
Gi . The M-sphere Bi needs not be real and may degenerate into a point; in the 
latter case, the cyclide is a ruled quadric (Bi = oo) or the image of a ruled 
quadric under an inversion. The M-sphere Bi may also be plane; in that case, 
the circle axes are the tangents of a conic in Bi or a pair of pencils of lines 
in Bi. 

The latter case appears when the hyperplane V* passes through the center Z 
of the stereographic projection and thus is mapped to a plane Bi in IR 3 . Now, 
Gi, the image of T* under the extended stereographic projection, is no longer 
a ruled quadric, but degenerates to the union of tangents of a conic or two 
pencils of lines. 

Thus we can construct paired circle families (and thus also Darboux cyclides) 
in a very simple way by prescribing a non-degenerate ruled quadric G and 
a sphere B (or sphere complex) and proceed according to Proposition M (see 



Fig. 11 and 12). Note that this construction would fail in the special case 
where B is a plane: prescribing a conic in a plane B is not enough to get the 
cyclide, since circles whose axes are tangent to the conic are automatically 
orthogonal to B. 
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Fig. 12. Various shapes of Darboux cyclides: A smooth Darboux cyclide is topolog- 
ically either a sphere (top), or a torus (middle), or a pair of two spheres (bottom 



right); see Takeuchi (2000). A cyclide with singularities in the image of a quadric 



under inversion (bottom left). Paired families of circles on the cyclides are shown. 

Remark 10 For readers familiar with non-Euclidean geometry we mention 
that Darboux cyclides may be obtained as follows: Take a ruled quadric in 
the projective model of hyperbolic or elliptic 3-space and transfer it (through 
the well-known Darboux transformation) to the respective conformal model. 
All cyclides carrying two paired families of circles which are orthogonal to a 
real or complex sphere can be obtained in this way; the former case is that of 
hyperbolic geometry, the latter belongs to elliptic geometry. 

Single circle families. If the cone Tj contains only one family of planes and 
thus has a vertex line its polar image is a conic T* =: q in the following 
sense: The points S* of the conic q determine the spheres S which envelope 
the cyclide D and the tangents of q represent the circles along which these 
spheres touch D. The plane P of q is the polar dual v* of Vf, it defines a bundle 
of spheres. Points on the vertex line Vi represent the spheres of a pencil, which 
are orthogonal to the spheres of the bundle. 

Recall that the spheres of a bundle which are orthogonal to an imaginary circle 
pass through two fixed points; this is the case of Darboux canal surfaces which 
are quadratic cones (if one of the two points is oo) or images of quadratic cones 
under an inversion. 

Proposition 11 A canal surface which is an irreducible Darboux cyclide (not 
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Fig. 13. A cyclide D which carries a single family of circles can be constructed as 
follows: Prescribe a conic g and a circle b in the same plane, and consider all spheres 
S which are centered on g and orthogonal to b. The cyclide is the envelope of these 
spheres. Each sphere S touches the cyclide in a circle k (of the single family), whose 
axis is a tangent of g and whose plane passes through the axis of b. 



of revolution) is the envelope of spheres S which are centered at a conic g and 
are orthogonal to an M-circle b in the plane of g (Fig. 13). The M-circle b 
may degenerate to a point and may also be imaginary; in the latter case, the 
cyclide is a quadratic cylinder (b — oo), a quadratic cone or an image of a 
quadratic cone or cylinder under an inversion. 



Remark 12 An exceptional case happens if the spheres are all centered on 
a straight line g (i.e., if the plane of the conic q passes through Z), i.e., the 
resulting surface is a rotational surface with axis g. The surface may be defined 
by rotating a circle around the axis, but only if the surface has another real 
circle family. There are cyclides for which this is not the case: as an example, 
take a rotational ellipsoid or the image of it under an inversion with respect 
to a co-axial sphere. We are not pursuing the study of cyclides with only one 
family of real circles, since we are interested in surfaces which carry at least 
two families of circles. 



3 Design tools for Darboux cyclides 



Section 2J3 shows that we can translate design methods for ruled quadrics and 
conies to design tools for paired and single circle families on Darboux cyclides. 
This will be elaborated in Section 3.1 in more detail. In Section 3.2 we give an 



example of circle families intersecting under constant angle and in Section 3.3 
we very briefly discuss the parametric representation of Darboux cyclides. 
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3.1 Designing circle families on cyclides. 



Most of the following material is based on Proposition [9j which generates 
paired families with help of a sphere complex and a ruled quadric. We will 
now illustrate this at hand of the case where the sphere complex consists of all 
spheres orthogonal to a real sphere B. It is pretty straightforward to discuss 
the other cases. 

The following design (input) methods for ruled quadrics will be transferred to 
cyclide design based on Proposition [9] 

(1) A ruled quadric is defined by three pairwise skew straight lines. 

(2) A ruled quadric is determined by a skew quad (formed by rulings) and a 
further point. 

(1) Translating the generation of a ruled quadric G with 3 pairwise skew 
rulings to cyclide design is done as follows: We prescribe a sphere B and three 
circles ki,k 2 ,k 3 which are orthogonal to B. Since the rulings of G shall be 
the axes of the circles, we have to prescribe the circles ki so that their axes 
Ai are pairwise skew. Two circles which are orthogonal to B have skew axes 
if and only if they do not lie on a common sphere. Having fixed the circles 
ki, the cyclide follows easily. One just has to compute a ruled quadric from 3 
pairwise skew lines Ai, A 2 , A 3 . One way is to first compute the rulings of the 
other family, as the set of lines intersecting A±, A 2 , A 3 . To get such a ruling, 
pick a point X G A\ and intersect the connecting planes X V A 2 and X V A$. 

It may be more intuitive to let B be a plane and let the user prescribe (half) 
circles orthogonal to B. Any two such circles should not be co-spherical. Lack- 
ing a direct construction in the case where B is a plane, one can first apply an 
inversion to map B to a sphere, proceed as above, and finally transform back 
with the (same) inversion. We will obtain a cyclide passing through the three 



prescribed circles (see Fig. 14). 



(2) The second type of input is closely related to the first one, but may be 
more useful for practice, since we now prescribe the circular boundaries and 



a point of the cyclide patch (see Fig. 15). Let us describe the case with a 



plane B. To this end, we can input a quad with vertices P\,...,P± and a 
plane B. Each edge, say P\,P 2 of the quad can be reflected at B to obtain 
points P{, P 2 , and now one has a unique circle passing through Pi, P 2 , P[, P 2 \ 
this circle ki is orthogonal to 5. So we end up with four circular arcs, being 
admissible boundaries of a cyclide patch. We may then prescribe one more 
point P of the desired patch and the cyclide D is determined. This input can 
be transformed into the above one by constructing a further circle k through 
P: take two opposite circular boundaries, say ki,k 3 and intersect the two 
connecting spheres P V ki and P V k 3 . 
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Fig. 14. Paired families of circles on a Darboux cyclide can be denned by 3 circles 
ki,k2, &3 which are orthogonal to a plane B (and no two of which lie on a common 
sphere). The figure shows a variety of results obtained with this method. 





Fig. 15. A cyclide patch can be defined by its four boundary circles (all being 
orthogonal to a common sphere or plane) and a further point. Especially if B is a 
horizontal plane, this can be very useful for architectural design, since one gets a 
structure formed by two families of circular arcs in vertical planes. 



Remark 13 Proposition 11 provides us with a transformation of conies into 
single families of circles on cyclides. Thus, in a quite analogous way to the 
design of paired circle families from ruled quadrics, we can transfer construc- 
tions of conies into those for single circle families on cyclides. This is rather 
straightforward, and thus we do not pursue it in the present paper. 



3.2 Two circle families intersecting under a constant angle 



Ivey (1995) proved that a surface which carries two orthogonal families of 



circles is a Dupin cyclide. Apart from the orthogonal network of principal 
curvature lines (circles) on a Dupin cyclide, there is also the case of the paired 
families of Villarceau circles on a ring cyclide (lying in doubly tangent planes) . 
Villarceu circles from different families always intersect under a constant angle, 
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as they intersect the principal curvature circles under a constant angle. For 
specific cyclides this angle can be 90 degrees (see Fig. 16). 



We see that any pair of circle families on a ring Dupin cyclide exhibits a 
constant intersection angle, which implies congruent nodes in an architectural 
construction. We conjecture that the Dupin cyclides are the only surfaces 
which have two circle families that intersect under a constant angle. 




Fig. 16. Since the Villarceau circles on a Dupin ring cyclide intersect the principal 
curvature circles under a constant angle, there are several ways of extracting circle 
families with a constant intersection angle. This leads to congruent nodes in an ar- 
chitectural construction. We show here two different types of 3-webs with congruent 



nodes, as discussed in Theorem 17 



3.3 Parametric representation of Darboux cyclides 



From the point of view of complex algebraic geometry a Darboux cyclide is 
a Del Pezzo surface of degree 4 (as a projection of the intersection of two 
quadrics in 4-space). In particular, Darboux cyclides are rational surfaces. 
However, in general they may not have a real rational parametrization. For 
instance, the cyclide (x 2 + y 2 + z 2 — l)(x 2 + y 2 + z 2 — 2) + 1/10 = is irreducible 
and disconnected, hence has no proper real rational parametrization. 

For an arbitrary Darboux cyclide D it is not hard to find a real parametrization 
of the form 



X(S t) = ^ +X2 ^ y[s t ) = Z{S t) ^ + ^ 



w l + w 2 vp w x + W 2 VP W 1 + w 2 vp 



where Xl, X 2 , Yi, Y 2 , Z±, Z 2 , W\, W 2 , P are polynomials in s and t of degree at 
most 4. 

Indeed, let X T -A-X = be the equation of D. Since the polynomial det(A-tJ) 
has degree 5, it has a real root tj. Thus D is the intersection of the cone 
X T • (A — tj J) • X = and the sphere E. Depending on the position of the 
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vertex V, of the cone with respect to the sphere £ there are the following 3 
cases. 

Case 1: Vi is in the "surface" of the sphere E. Consider the hyperplane Hi 
touching the sphere £ at the point opposite to Vi. The hyperplane Hi intersects 
the cone by a quadric Di. Take a quadratic parametrization of the quadric Di. 
Composing it with the inverse stereographic projection from the vertex Di 
and then with the stereographic projection from the point Z we get a rational 
parametrization of the cyclide D of degree a most 4. 

Case 2: Vi is outside the sphere S. Perform a projective transformation k of 
the 4-space preserving £ and taking V to the infinitely distant point in the 
direction of the xi-axis. The image of the cone under the transformation k, 
is a cylinder crossing the hyperplane X\ = by a quadric (or a plane) Di. 
Thus the quadric (or a plane) Di is the projection of kD along the xi-axis. 
Take a quadratic parametrization of the quadric Di. Inverting the projection, 
performing the inverse transformation and finally performing the stereo- 
graphic projection from Z we get the required parametrization. 

Case 3: Vi is inside the sphere £. Perform a projective transformation k of 
the 4-space preserving £ and taking V to the center of the sphere £. The 
image of the cone under the transformation intersects the hyperplane x\ = 1 
by a quadric (or a plane) Di. Thus the quadric (or a plane) Di is the central 
projection of the cyclide from the origin. Take a quadratic parametrization of 
the quadric Di. Inverting the projection, performing the inverse transformation 
k~ x , and finally performing the stereographic projection from Z we get the 
required parametrization. 



4 Webs from circles on Darboux cyclides 



In this section we define 3-webs, give classical examples of planar C-webs, and 
finally classify all possible C-webs on Darboux cyclides. 

A 3-web (or hexagonal web) in a surface is 3 families of smooth curves which 
are locally diffeomorphic to 3 families of lines x = const, y = const, x + y = 
const in the plane Oxy. These lines for integral values of the constants form 
a triangular lattice in the plane; taking the preimage of the lattice under the 



diffeomorphism, we get a discrete 3-web in the surface (see Fig. 18). 



A 3-web is characterized by the following closure condition (Blaschke and Bol 



(1938); see also Fig. 17). Take an arbitrary point O of the surface. Draw the 
curves k\, k 2 , k% of the 3 families through the point O. Take an arbitrary 
point A 1 G ki sufficiently close to O. Draw the curve l 2 of the second family 
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Fig. 17. The closure condition characterizing 3- webs; see the beginning of Section [4] 
for explanation. This figure also illustrates an algorithm for computation of discrete 
3- webs on Darboux cyclides; see the end of Section [4} 

through the point A\. Set A 2 = h H k%. Draw the curve Zi of the remaining 
(first) family through the point A 2 . Set A3 = l\ n k 2 . Draw the curve I3 of 
the remaining (third) family through A 3 and set A4 = l 3 D k±. Continuing in 
this way, construct points A 5 , A 6 , and A 7 . The closure condition asserts that 
A-j = A\. This condition and symmetry considerations often allow to verify 
that certain triple of families is indeed a 3-web. 

First let us give the following known beautiful examples of 3-webs of lines 
and circles in the plane needed for the sequel. In these examples by a 3-web 
we mean a 3-web in an appropriate nonempty domain in the plane. Three 



families of curves are 



crossing, 



if some 3 curves of distinct families have a 



common point of transversal intersection. 



Theorem 14 (Graf and Sauer (1924); Blaschke and Bol (1938)) Three cross- 
ing families of (pairwise distinct) straight lines in the plane form a 3-web if and 
only if all these lines are tangent to one curve of algebraic class 3 (including 
reducible ones). 

A cyclic is a curve given by equation of the form 

X(x 2 + y 2 ) 2 + (x 2 + y 2 )(fix + vy) + Q(x, y) = 0, 

where A,//, v G K and Q(x,y) is a polynomial of degree at most 2. (This is a 
one-dimensional counterpart of a cyclide.) 



Theorem 15 (Wunderlich (1938)) Three crossing families of circles doubly 
tangent (in the sense of complex algebraic geometry) to a cyclic (possibly with 
singularities) form a 3-web. 
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Remark 16 Performing the inverse stereographic projection and then a pro- 
jective transformation one gets the following corollary: the tangent planes to 3 
of the 4 quadratic cones (possibly reducible) passing through the intersection c 
of two oval quadrics in M. 3 cut a 3-web of conies in each of the oval quadrics, 
once the 3 families of conies are crossing (cf. Fig. [?[). 

Since we have up to 6 real families of circles on a cyclide, the question arises 
whether and how we can arrange triples of them into 3-webs. To this end, we 



will prove the following result (see Fig. 18): 



Theorem 17 Three families of circles on a nontrivial irreducible Darboux cy- 
clide form a 3-web unless one takes two nonspecial paired families and another 
family which has a paired one. Thus we have 5 types of 3-webs from circles on 
a Darboux cyclide: 



(i) A web of type 1 is formed by 3 non-single families of circles, such that no 
two of them are paired families. Hence, this type only exists on a cyclide 
with 6 real families of circles, and we have 8 different 3-webs on such a 



cyclide; see Fig. 18 



(ii) A web of type 2 is formed by 2 special paired families and another family 
which has a paired one. This type only exists on the image of a one-sheet 
hyperboloid under an inversion; see Fig. [Tg[ 

(Hi) A web of type 3 is formed by a single family and 2 paired families; see 
Fig. \fi^ to the left. 

(iv) A web of type 4 is formed by a single family and 2 non-single families, 
which are not two paired families; see Fig. 

(v) A web of type 5 consists of 2 single families and another family which has 
a paired one. Hence, this type only exists on the Dupin ring cyclide and 
consists of its principal curvature lines and a family of Villarceau circles; 



see Fig. 16 to the right. 



Proof. The main idea of the proof is to use the surface L) C M 4 , the cones 
which define the circle families and to verify webs at hand of appropriate 
projections onto a plane or 3-space. Types 2,4,5 are in fact degenerations of 
type 1. 



Type 1: Here D is the carrier of a pencil of quadrics which contains three 
cones (with vertices Vi,V2,V^) each of which contains planes. By Remark [8] 
the pencil of quadrics must also contain a cone with a real vertex V4 and an oval 
base quadric. Now we project D from V4 onto a hyperplane if 4 (not through 
V4) and obtain (a part of) an oval quadric $. The three circle families of L) 
(no two of them paired) are projected to conies (on $) in tangent planes of 
quadratic cones; $ and these cones lie in a pencil of quadrics. By Proposition [5] 
it follows that the conies of distinct families are crossing. We have now exactly 
the situation of Remark [16] and conclude that these three families of conies 
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18. Eight different 3-webs of type 1 in a Darboux cyclide. 
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Fig. 19. A 3-web of type 2 in a Darboux cyclide (right) obtained from a non-rota- 
tional one-sheet hyperboloid (left) by an inversion. 

form a 3-web. Hence, also the 3 circle families on the cyclide form a 3-web. 

Type 2: Let V\ G £ be vertex of the cone giving the special paired families. Let 
V2 be the vertex of the cone Y2 giving the other non-single family. Project from 
the line connecting V\ and V2 onto a plane. Then the special paired families 
are projected either to tangents to a conic or to two pencils of lines. Since 
V\ G £ it follows that V\ G T 2 . Thus the 3rd non-single family is projected to 
a pencil of lines. The resulting straight lines are tangent to a reducible curve 



of algebraic class 3 and by Theorem 14 form a 3-web. Hence, also the 3 circle 
families on the cyclide form a 3-web. 

Special paired families only exist on cyclides which are images of quadrics 
under inversion. Among the quadrics only a one-sheet hyperboloid (not of 
revolution) contains 4 families of circles/lines. Thus the cyclide in type 2 can 
only be an image of a one-sheet hyperboloid under inversion. (This web on a 



hyperboloid was also discussed by Wunderlich (1938).) 



Type 3: We connect the cone vertex V\ (associated with the two paired families) 
with any point on the vertex line v 2 of the cone that determines the single 
family. From the resulting line L we project into a plane (not intersecting L). 
The image of the two paired circle families are either the tangents of a conic 
or two pencils of lines. The image of the single family is a pencil of lines. So 
together we have the tangents of a reducible curve of algebraic class 3, which 



according to Theorem 14 form a 3-web 



Type 4: The proof is literally the same as in the proof for type 1. We only need 
to remark that one of the cones in the hyperplane H4 now degenerates to a 
pair of planes and its tangent planes are understood in the sense of algebraic 
geometry: these are the ones passing through the vertex line of the cone. 

Type 5: We may apply an inversion to map the ring Dupin cyclide to a ring 
torus. Now it follows from the rotational symmetry that the Villarceau circles 
together with meridian circles and parallel circles (rotational paths) form a 
3-web. 
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Fig. 20. A 3-web of type 4 in a Darboux cyclide. 



We still have to show that two (nonspecial) paired families and a 3rd (non- 
single) circle family cannot form a 3-web. Again, project from the line con- 
necting the two involved cone vertices, say Vi, V2 onto a plane. The two paired 
families are projected either to tangents to a conic or to two pencils of lines. 
Since V\ S it follows that the 3rd family is projected to tangents to a (non- 
degenerate) conic. The resulting straight lines are tangent to a reducible curve 



of algebraic class 4 and due to Theorem 14 do not form a 3-web. □ 



Remark 18 Similarly, there are the same 5 types of 3-webs of conies in the 
intersection of two arbitrary 3 -dimensional nonsingular quadrics in 1R 4 (if the 
intersection does not contain 2-dimensional quadrics as a whole). If we project 
such a 3-web from M 4 to IR 2 we get a 3-web of conies in the plane. The conies 
of the planar 3-web are tangent to the "visual contour" (i.e., the boundary of 
the projection) of the intersection of the quadrics, which is a curve of order 



< 12; cf. Schmid (1936), Timorin (2007). 



The closure condition can be used for the computation of discrete 3-webs as 
follows. Consider 3 families of curves on a surface, forming a 3-web. We use 
the notation from the definition of the closure condition at the beginning of 
Section [4} We prescribe the requested "size" N of the web, as well as starting 
points O and Ai 6 k%. One step of algorithm is to draw the curve m 3 of the 
third family through the point Ai, set A\ = m 3 R l\\ then draw the curve m 2 
of the second family through the point A' x \ set A' 2 = PI k%. The next step 
of algorithm is the same, but now A2 and A\ are chosen as starting points. 
Performing N such steps one gets A" curves of the first family and A" curves of 
the second one. After that, we make similar A" steps starting from the points 
O and A 2 to get A" curves of the third family. 

One can apply this algorithm to construct finite 3-webs, i.e., discrete 3-webs 
forming a triangulation of the whole closed surface. Indeed, parameterize one 
of the curve families, say, the first one, so that the curve k\ corresponds to 
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parameter values both and 1 (we assume that the surface has ring shape). 
Running the algorithm we get N curves of the first family, let a be the param- 
eter value of the iV-th curve. Let A\ e k\ be a point such that OA\ = OA\/o~. 
Repeat the same procedure starting from the points O and A\ until a is close 
enough to 1. Eventually a discrete 3- web close to a finite one is constructed 
this way. On Darboux cyclides, there are arbitrary dense finite 3-webs of type 
1 by dWunderlichj |1938[ Section 4). 



5 Related studies and open problems 



There are several directions in which we would like to continue our research: 

• We are still lacking a proof for the conjecture that any surface which carries 
three families of circles is a Darboux cyclide and that any surface with two 
families of circles intersecting under a constant angle is a Dupin cyclide. 
Moreover, we would like to classify all surfaces which carry two families of 
circles and discuss them from a geometric modeling perspective. 

• Cyclides of isotropic geometry may also deserve attention. Recall that an 
isotropic circle is either a parabola with a vertical axis (a favorite element 
of architecture) or an ellipse whose projection into a horizontal plane (top 
view) is a circle. An isotropic cyclide is given by equation ([I]), in which 
both instances of (x 2 + y 2 + z 2 ) are replaced by (x 2 + y 2 ). Viewing isotropic 
cyclides from the Laguerre geometric perspective, we can transform them 
into surfaces which are enveloped by several families of right circular cones. 
This would, for example, provide access to architectural structures from 
strips of right circular cones. 

• It is also interesting to obtain an algorithm for rational parametrizations of 
Darboux cyclides and the corresponding rational Bezier representation; this 



is expected to result from the work of Krasauskas and Zube (2011). 
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